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Data  indicate  that  certain  important  types  of  marine  organisms  behave  acoustically  like  weakly 
scattering  fluid  bodies  (i.e.,  their  material  properties  appear  fluidhke  and  similar  to  those  of  the 
surrounding  fluid  medium).  Use  of  this  boundary  condition,  along  with  certain  assumptions, 
allows  reduction  of  what  is  a  very  complex  scattering  problem  to  a  relatively  simple, 
approximate  ray-based  solution.  Because  of  the  diversity  of  this  problem,  the  formulation  is 
presented  in  two  articles:  this  first  one  in  which  the  basic  physics  of  the  scattering  process  is 
described  where  the  incident  sound  wave  is  nearly  normally  incident  upon  a  single  target  (i.e., 
the  region  in  which  the  scattering  amplitude  is  typically  at  or  near  a  maximum  value  for  the 
individual)  and  the  second  one  [Stanton  et  al.,  J.  Acoust.  Soc.  Am.  94,  3463-3472  (1993)] 
where  the  formulation  is  heuristically  extended  to  all  angles  of  incidence  and  then  statistically 
averaged  over  a  range  of  angles  and  target  sizes  to  produce  a  ccllcctive  echo  involving  an 
aggregation  of  randomly  oriented  different  sized  scatterers.  In  this  article,  a  simple  ray  model  is 
employed  in  the  deformed  cylinder  formulation  [Stanton,  J.  Acoust.  Soc.  Am.  86,  691-705 
( 1989)]  to  describe  the  scattering  by  finite  length  deformed  fluid  bodies  in  the  genera)  shape  of 
elongated  organisms.  The  work  involves  single  realizations  of  the  1  ngth  and  angle  of 
orientation.  Straight  and  bent  finite  cylinders  and  prolate  spheroids  are  treated  in  separate 
examples.  There  is  reasonable  qualitative  comparison  between  the  structure  of  the  data  collected 
by  Chu  et  al.  [ICES  J.  Mar.  Sci.  49,  97-106  (1992)]  involving  two  decapod  shrimp  and  this 
single-target  normal-incidence  theory.  This  analysis  forms  the  basis  for  successful  comparison 
(presented  in  the  companion  article)  between  the  extended  formulation  that  is  averaged  over  an 
ensemble  of  realizations  of  length  and  angle  of  orientation  and  scattering  data  involving 
aggregations  of  up  to  lOO’s  of  animals. 

PACS  numbers:  43.20.Fn,  43.30.Gv,  43.30.Xm 
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scattered  pressure 

plane-wave/plane  interface  reflection  coefficient 
where  wave  is  initially  in  medium  i  and  in¬ 
cident  upon  medium  j  [=(pjCj/p/Ci—l)/ 
iPjC/PEi+^)] 

transmission  coefficient  for  planar  interface  due  to 
plane  wave  initially  in  medium  i  and  in¬ 
cident  upon  medium  j  [Tij=2(pjCj/p/Ci)/ 
(l  +  (pf/PiCi))] 

subscripts  indicating  medium  “1”  (surrounding 

fluid)  and  medium  “2”  (body  medium) 

acoustic  wave  number  ( =2i7/A) 

radius  of  cylinder  cross  section 

amplitude  of  incident  plane  wave 

coefficient  determined  from  boundary  conditions 


4>  azimuthal  angle  {d>  =  iT  is  the  backscatter  dire^T^ 
tion)  t 

/ 

form  function  for  infinitely  long  cylinder 
/(k^a)  scattering  amplitude  for  finite-sized  objects  (sub¬ 
scripts  sc,  ps,  and  be  refer  to  straight  cylinder, V 
prolate  spheroid,  and  bent  cylinder) 
r  distance  from  axis  of  cylinder  to  the  field  point 

(applied  to  infinite  and  finite  bodies) 
fpos  position  vector  from  the  origin  of  the  coordinate  ^ 
system  to  a  point  on  the  axis  of  the  cylinder  Q 
r,  distance  from  a  point  on  the  axis  (at  rp„J  to  the 
field  point  (receiver) 

ap  term  in  divergence  factor  for  scattered  ray 

Bp  reflection-transmission  factor 

T]p  propagation  pliase  delay  of  scattered  ray 
Pp  phase  advance  associated  with  crossing  of  caustics  ^ 
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‘  fig  phase  shift  associated  with  external  caustics 

tj  distance  between  the  point  on  the  axis  (at  fpos) 

and  the  plane  that  both  contains  the  origin  and  is 
perpendicular  to  the  direction  of  the  incident 
plane  wave 
A  A:,Z,sin^ 

q  apparent  volume  flow  per  unit  length  of  scattered 

field  due  to  infinitely  long  cylinder 
6  angle  between  direction  of  incident  plane  wave 

and  plane  whose  normal  is  the  tangent  to  the  axis 
at  each  point  (“angle  of  incidence”  to  local 
tangent) 

f,,f,  unit  vectors  describing  direction  of  incident  and 
received  (scattered)  waves 

L  length  of  straight  finite  object;  arc  length  (or  pro¬ 
jected  length,  depending  on  context )  of  uniformly 
bent  cylinder 

Lebc  effective  length  of  bent  cylinder 


INTRODUCTION 

Quantitative  remote  sensing  of  marine  organisms  with 
sonars  requires  detailed  knowledge  of  their  scattering  prop¬ 
erties.  In  a  series  of  papers,  the  deformed  cylinder  formu¬ 
lation  has  demonstrated  promise  for  the  description  of  the 
scattering  by  certain  elongated  zooplankton  and 
micronekton'"*  and  fish.’  In  particular,  straight,  bent,  and 
rough  cylinder  evaluations  of  the  formulation  have  de¬ 
scribed  some  or  most  of  the  scattering  characteristics  of 
shrimplike  (crustacean)  zooplankton — an  abundant  class 
of  animals.  The  analyses  suggest  that  most  of  these  ani¬ 
mals,  whose  bodies  tend  to  be  bent,  behave  acoustically  like 
bent  cylinders’"*  although  some  data  imply  that  certain 
animals  may  behave  like  straight  cylinders.'  The  rough 
cylinder  formulation®  ’  has  described  the  statistical  nature 
of  the  scattering  by  moving,  flexing  animals.®  Clay  has 
successfully  described  the  scattering  of  sound  by  fish  by 
modeling  the  swimbladders  as  bent  gas-filled  cylinders  and 
the  fleshy  part  of  the  fish  as  weakly  scattering  bent  fluid 
cylinders.’  In  that  paper,  he  showed  that  the  apparent 
“damping”  effect  at  resonance  hypothesized  by  others  can 
be  explained  by  the  elongated  shape  of  the  nonspherical 
swimbladder. 

To  date,  most  of  the  uses  of  the  deformed  cylinder 
formulation  have  involved  the  modal  series  solution  which 
requires  many  terms  to  converge  in  the  high-frequency  re¬ 
gion  and,  in  general,  is  very  cumbersome,  if  not  impossible, 
to  manipulate  algebraically.  Manipulation  is  essential  when 
deriving  formulas  (usually  approximate)  that  not  only  ex¬ 
plicitly  illustrate  the  physics  of  the  scattering  process  but 
also  describe  the  average  scattering  properties.  For  exam¬ 
ple,  in  order  to  estimate  the  scattering  by  rough  cylinders, 
an  approximate  ray  solution  was  used  in  the  formulation.®*’ 
The  results  were  used  to  describe  the  general  trend  of  the 
statistics  of  the  echo  from  zooplankton®  although  since 
that  solution  involved  dense  elastic  objects,  it  was  not  ap¬ 
propriate  to  compare  the  predicted  levels  to  those  mea¬ 
sured  from  the  animals.  In  addition  to  the  rough  cylinder 
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Oq  length  of  semi-minor  axis  of  prolate  spheroid 
Pg  radius  of  curvature  of  axis  of  uniformly  bent  cyl¬ 
inder 

p  mass  density 

Y  position  angle  of  bent  cylinder 

Tniax  angle  that  subtends  portion  of  bent  cylinder  be¬ 
tween  midpoint  and  end 

X  position  of  bent  cylinder  axis  that  has  been  pro¬ 
jected  onto  X  axis 

A  (maximum  deflection  of  bent  cylinder) /a 

/?  (length  of  cylinder)/a 

^  Fresnel  integral  parameter 

C,5  Fresnel  integrals 

c  compressional  speed  of  sound 

8  Pi/Px 

h  c-i/c^ 

RTS  target  strength  —10  log/,’  (reduced  target 
strength) 


problem,  simple  solutions  are  necessary  to  describe  average 
echo  levels  from  aggregations  of  animals  of  random  orien¬ 
tation  and  length. 

Data  collected  over  the  past  several  years  provide  in¬ 
formation  allowing  us  to  develop  a  simple  model.  In  par¬ 
ticular,  field  data  collected  by  Pieper  et  aL,^°  and  labora¬ 
tory  data  collected  by  Stanton  eta/."  involving  average 
echoes  from  aggregations  show  a  dip  in  the  backscatter 
versus  ka  curve  toward  the  lower  end  of  the  geometric 
scattering  region.  Laboratory  data  from  Chu  et  al}  involv¬ 
ing  echoes  from  each  of  two  individuals  show  a  deep  null  in 
the  same  region.  The  laboratory  data  presented  in  the  latter 
two  articles  showed  the  dip  or  null  to  consistently  lie  in  the 
region  surrounding  ka=2  where  a  is  the  equivalent  cylin¬ 
drical  radius  of  the  elongated  animals  (decapod  shrimp 
whose  bodies  are  representative  of  a  broad  class  of  marine 
organisms).  The  dip  or  null  is  indicative  of  strong  inter¬ 
ference  between  the  acoustic  wave  or  “ray”  reflected  from 
the  front  interface  of  the  animal  and  a  ray  that  has  pene¬ 
trated  the  interior  with  little  loss  and  reflected  off"  the  back 
interface.  Complementing  these  data  are  direct  measure¬ 
ments  of  density  and  sound  speed  of  various  zooplankton 
that  indicate  that  those  properties  are  very  similar  (to 
within  several  percent)  to  those  of  the  surrounding 

12-14 

seawater. 

All  of  the  above  (direct  measurement  and  scattering) 
data  viewed  collectively  suggest  that  the  animals  tend  to 
behave  acoustically  as  weakly  scattering  bodies.  For  a  body 
to  be  weakly  scattering,  the  material  properties  must  be 
similar  to  those  of  the  surrounding  medium  which,  as 
stated  above,  is  the  case  for  these  zooplankton.  Also,  as 
will  be  illustrated  in  the  theoretical  and  numerical  portions 
of  this  article,  weakly  scattering  bodies  of  certain  classes  of 
shapes  (cylindrical  in  this  case)  can  be  described  acousti¬ 
cally  by  a  simple  two-ray  model  at  broadside  incidence. 
Such  a  formulation  gives  rise  to  a  deep  null  near  ka=2  for 
a  single  target  which  is  consistent  with  the  data  collected 
by  Chu  eta/.®  Note  that  even  the  average  echo  from  ag¬ 
gregations  of  randomly  oriented  animals  contains  a  dip 
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that  appears  to  be  a  “smeared”  version  of  the  ‘  The 

dip  is  strong  enough  (of  the  order  5-10  dB  below  sur¬ 
rounding  levels)  that  a  model  needs  to  be  developed  that 
predicts  it. 

With  this  significant  evidence  that  the  animals  behave 
as  weakly  scattering  bodies,  we  formulaie  a  simple  analyt¬ 
ical  approximate  model  that  estimates  the  scattering  levels 
and  illustrates  the  fundamentals  of  the  scattering  process. 
Choice  of  this  model  was  based,  in  part,  on  the  desire  for 
one  that  can  be  easily  manipulated  algebraically  and  is 
computationally  efficient.  These  last  three  desirable  fea¬ 
tures  are  not  possible  with  the  modal-series-based  model 
described  in  Ref.  2  or  other  numerically  oriented  models 
although  note  that,  while  this  current  model  is  limited  to 
weakly  scattering  materials,  the  modal-series-based  model 
can  be  used  for  a  broad  range  of  materials. 

The  analysis  is  divided  into  two  articles:  this  present 
one  which  involves  formulating  the  basic  physics  of  the 
scattering  process  for  a  single  realization  of  length  and 
angle  of  orientation  in  the  region  where  the  scattering  am¬ 
plitude  tends  to  be  at  a  maximum  (near  normal  incidence) 
and  a  companion  paper  where  the  work  is  heuristically 
extended  to  include  all  angles  of  incidence  for  the  purpose 
of  averaging  over  angle.**  The  resultant  extended  model  is 
averaged  in  Ref.  11  over  a  range  of  angles  and  sizes  to 
produce  a  model  describing  an  “aggregation  echo.”  The 
aggregation  model  is  successfully  compared  in  Ref.  1 1  with 
data  involving  aggregations  of  animals  ranging  in  numbers 
from  six  to  many  hundreds.  The  data  show  the  same  5-  to 
10-dB  dip  as  predicted  by  the  model,  hence  further  vali¬ 
dating  the  ray-based  approach  described  in  this  present 
paper. 

In  this  present  article,  a  simple  existing  ray  model  is 
modified  and  incorporated  in  the  deformed  cylinder 
formulation.^  The  modification  involves  heuristically  alter¬ 
ing  the  phase  of  the  rays  so  that  the  ray  description,  which 
is  normally  valid  only  in  the  geometric  or  ka>  1  region,  can 
be  used  accurately  for  values  of  ka  well  below  unity.  While 
the  derivation  is  motivated  by  the  animal  problem,  the 
results  are  general  and  not  specific  to  scattering  by  animals. 
Examples  involving  several  shapes,  comparison  with  the 
modal-series-based  solution  that  involves  fewer  mathemat¬ 
ical  approximations,  and  comparison  with  data  from  single 
animals  are  given. 

I.  THEORY 

A.  Oerwral  ray-bMed  solution 

The  far-field  solution  to  the  scattering  by  an  infinitely 
long  cylinder  is  given  as*^ 

kr>l,  (1) 

where  the  form  function  is 
2  ^ 

X  (2) 

yjirka  m-o 

The  modal  series  representation  in  Eq.  (2)  is  exact  for 
all  ka  and  material  composition  (axisymmetric  profile). 
For  kayl,  however,  many  terms  are  required  for  mathe- 

3466  J.  Acouat  Soc.  Am.,  Vd.  94,  No.  6,  December  1993 


Incidsnt  pione  wavs 


no.  1.  (t)  Weak  scattcrer  geometry  showing  two  types  of  rays  scattered: 
“specular"  echo  (/,.o)  reflected  olT  front  interface  and  "transmitted” 
echo  reflect^  ofl'back  interface  after  having  traveled  through  the 

body,  (b)  Deformed  cylinder  and  general  bistatic  sonar. 

matical  convergence  and,  while  modem  computers  can  ef¬ 
ficiently  calculate  the  series  to  high  precision,  the  series  is 
quite  cumbersome  to  manipulate  algebraically.  Use  of  ray 
solutions  in  this  geometric  scattering  region  has  proven  to 
facilitate  the  manipulations  while  at  the  same  time  making 
the  physics  of  the  scattering  process  more  explicit  in  the 
mathematics.’"’’ **“^*  The  ray  theories  are  based  on  a  vari¬ 
ety  of  methods  including  the  Sommerfeld-Watson  trans¬ 
formation,  heuristically  derived  formulations,  and  the 
Kirchhoff  method.  The  success  of  using  ray  methods  to 
describe  the  scattering  by  marine  organisms  is  illustrated 
by  Foote^*  where  the  Kirchhoff  method  was  used  to  nu¬ 
merically  calculate  the  scattering  by  the  front  interface 
alone  from  (gas-filled)  swimbladders  of  fish. 

The  types  of  rays  that  need  to  be  taken  into  account 
depend  upon  the  material  properties.  In  the  weak  fluid 
scatterer  case  analyzed  in  this  article  (g,h~  1),  it  is  shown 
that  only  two  rays — the  ones  reflected  from  the  front  and 
back  interfaces — need  to  be  taken  into  account  for  reason¬ 
able  (although  not  exact)  description  of  the  scattering  in 
the  transitional  and  geometric  region  [Fig.  1  (a)].  The  form 
function  associated  with  each  ray  can  be  calculated  from 
the  general  (geometric  scattering)  formula  for  the pth  scat¬ 
tered  ray  from  Marston:’* 

=  [2/u ) I  oip-a\  (3) 
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where  thep=0  and  />=2  rays  correspond  to  the  reflections 
off  the  front  and  (first-order  reflection)  back  interfaces, 
respectively.  The  values  oi  p>  2  correspond  to  higher- 
order  internal  reflections.  According  to  the  formulas  in 
Marston,  ap._^Q=a/2,  —Ik^a,  and 

for  the  p  =  0  “specular”  reflection  off  the  front 
interface  while  0^=2=  (3/2 )a  (weak  scatterer  approxima¬ 
tion),  B=  —  7’,27'2,."s/i2,  7jp=2  =  2k,a(2A~'- 1),  and 

/i^^2=— for  the  p  =  2  “transmitted  wave”  reflection 
off  the  back  interface.  Here,  /i,=0  for  all  rays  in  the  case 
presented  in  this  article  where  C2>Ci. 

Although  Eq.  (3)  is  generally  valid  only  in  the  geo¬ 
metric  optics  limit,  we  have  found  that  its  usefulness  can  be 
extended  down  to  the  Rayleigh/geometric  transitional  re¬ 
gion  (0. 1  S  /ca  1 )  by  heuristically  removing  the  phase  ad¬ 
vance  fXp^2  that  is  due  to  a  caustic.  This  caustic  is  related 
to  the  curvature  of  the  cross  section  of  the  body.  Phase 
shifts  due  to  the  caustic  naturally  reduce  to  zero  in  the  long 
wavelength  limit.  While  the  variation  of  phase  with  respect 
to  wavelength  can  be  derived,  we  remove  the  caustic  effect 
gradually  by  replacing  the  constant  value  —'^/2  by 

the  convenient  function  Pp^iikxa)  =  ~{'n/2)kxa/ 
(/cia-)-0.4).  The  high-frequency  (A:ia>l)  and  low  fre¬ 
quency  (^ia<l)  limits  of  Pp^iik^a)  are  —ir/2  and  0,  re¬ 
spectively.  The  constant  0.4  was  determined  empirically 
and  limits  the  errors  to  within  about  2  dB  for  the  values 
A:u>0.l  and  (g,/i)<(l.l,l.l)  for  the  straight  cylinder. 

Evaluating  Eq.  (3)  with  these  parameters  gives  the 
approximate  form  function  for  the  weakly  scattering  infi¬ 
nitely  long  fluid  cylinder: 

(4) 

=^,2e"^*'‘'(  I  -  r,2r2iC^*^V^''=2^*'‘"),  (5) 

where  the  general  form  for  the  scattering  is  given  in  Eq. 
(4)  and  an  explicit  approximate  form  is  in  Eq.  (5).  The 
assumptions  for  these  equations  are  ( I )  the  incident  plane 
wave  is  traveling  in  a  direction  normal  to  the  axis  of  the 
cylinder,  (2)  the  cylinder  is  surrounded  by  a  fluid,  and  (3) 
the  scattering  is  sufficiently  weak  that  (a)  the  first  reflec¬ 
tion  off  the  back  interface  dominates  the  higher-order  in¬ 
ternal  reflections  and  (b)  bending  of  the  rays  at  the  inter¬ 
face  (Snell’s  law)  is  negligible.  Simulations  later  in  this 
article  verify  the  validity  of  assumption  (3). 

A  useful  property  of  the  scattered  field  from  elongated, 
deformed,  sometimes  finite  bodies  is  the  apparent  volume 
flow  per  unit  length  q  of  the  field. This  term  can  be 
integrated  over  the  length  of  the  body  to  provide  estimates 
of  the  scattered  field,  and  hence  backscattering  cross  sec¬ 
tion  and  target  strength.  The  accuracy  of  the  estimates 
improve  as  the  aspect  ratio  (length/diameter)  increases. 
The  volume  flow  per  unit  length  will  be  used  in  predicting 
scattering  by  the  finite  objects  in  this  article. 

The  far-field-scattered  pressure  due  to  a  plane-wave 
incident  upon  an  infinitely  long  cylinder  is  written  in  terms 
of  9  as 

srr)qe-‘'"*  \'^p,c',/(2  ^/^)  (/cir>  I), 

(6) 


where 


_  4Po 


X  bm  COS  m({>. 


m~0 


(7) 


The  volume  flow  per  unit  length  can  be  related  to  the 
form  function  with  the  following  expression: 


q=2Po^{  -rra/k  “  VpiC, .  (8) 

Inserting  this  equation  in  Eqs.  (4)  and  (5)  gives  an  ap¬ 
proximate  expression  for  the  volume  flow  per  unit  length 
due  to  the  scattering  from  a  weakly  scattering  fluid  cylin¬ 
der: 


Q  —  <lp=0  +  9p=lt 


(9) 


=  2Fo^i2c"^*'‘'(  I  -  7’,27’2ie^*2V>'^=2(*i“)) 


Vto/(piC,  ^).  (lO) 

The  above  volume  flow  term  is  incorporated  into  the 
deformed  cylinder  formulation  so  that  effects  due  to  defor¬ 
mities  such  as  bend,  taper,  and  roughness  can  be  estimated. 
The  approximate  formulation  allows  for  calculation  of  fi¬ 
nite  or  infinitely  long  bodies.  The  general  solution  for  bod¬ 
ies  of  any  length  is  given  as 


Pscal 


q  cos  6  ■ 


I  dr, 


po$l 


X(any  length),  (ll) 

where  some  of  the  terms  are  illustrated  in  Fig.  l(b).  All 
wavenumbers  in  q  are  modified  by  the  factor  cos  6.  Note 
that  in  Eq.  (4)  of  Ref.  2  there  is  no  explicit  cos  9  term  that 
corresponds  to  the  one  that  appears  in  the  above  expres¬ 
sion.  The  misleading  omission  in  Ref.  2  made  no  difference 
in  the  results  because  from  Eq.  (S)  and  beyond,  the  modal 
series  representation  of  q  was  used  and  a  ( cos  0)  ” '  term  in 
q  cancelled  out  the  would-be  cos  6  factor.  In  the  ray  rep¬ 
resentation,  there  is  no  such  cancellation  hence  the  cos  6  is 
written  explicitly.  The  solution  for  bodies  of  finite  length 
(i.e.,  the  length  is  much  less  than  the  first  Fresnel  zone  of 
the  source/receiver  combination)  simplifies  to 

Psca,='Po(^'*'Vr)/(k,o),  (I2a) 

where 


/(M)  =  (^o) 


q  cos  6 


Xexp[/k,rpo,(f,-r,)  •  Tpos]  |drpos| 

(finite  length).  (I2b) 

The  details  of  these  equations  are  discussed  in  Refs.  2 
and  8.  Note  that  the  implicit  Pq  in  q  cancels  the  explicit 
(Po)“‘ inEq.  (I2b). 


B.  Deformed  smooth  bodies 

The  ray  formulation  can  readily  be  incorporated  into 
the  above  deformed  cylinder  equations  by  direct  substitu¬ 
tion  to  obtain  approximate  solutions  to  the  scattering  by  a 
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variety  of  bodies  such  as  the  straight  finite  cylinder,  prolate 
spheroid  (high  aspect  ratio),  and  uniformly  bent  cylinder. 
Adaptation  of  the  integral  formulation  to  each  shape  is 
described  in  Refs.  2  and  8  and  includes  ( 1 )  allowing  the 
“cylindrical  radius”  of  the  prolate  spheroid  to  vary  along 
the  lengthwise  axis,  (2)  taking  into  account  the  phase 
shifts  induced  by  the  bend  of  the  axis  of  the  bent  cylinder 
(this  takes  place  in  the  exponent  in  the  integrand),  and  (3) 
solving  the  integral  along  the  length  of  the  cylinders  for  the 
prolate  spheroid  and  bent  cylinder  by  use  of  the  method  of 
stationary  phase  or  Fresnel  integrals.  Performing  the  above 
manipulations  gives  the  following  approximate  solutions 
for  those  shapes  in  the  geometric  scattering  region; 

1.  Straight  cylinder  (various  angles  of  incidence); 

X(l-r,ir2,e'^2'‘=“*e'^*>=2'*'‘'').  (13) 

2  Prolate  spheroid  (broadside  incidence): 

/ps“  ( 1/4)  ue-'’*'‘^(  ). 

(14) 

3.  Bent  cylinder  (uniformly  bent  or  nearly  so;  bent 
symmetrically  away  from  sonar). 

a.  Arbitrary  deflections.  For  values  of  deflection  of  the 
bent  cylinder  that  are  arbitrary  with  respect  to  wavelengths 
of  sound,  one  obtains 

“  ( 15a) 

»=2ebc/«(^=0*)/2,  (15b) 

where  the  approximations  ki^ki  and  cos  I  were  made 
in  amplitude  terms  and  the  limit  Pc>o  was  used  (which  is 
reasonable  for  elongated  zooplankton).  The  effective 
length  can  be  expressed  in  terms  of  either  the  radius  of 
curvature  of  the  cylinder  axis  or  deflection 
z(x)  =4Ax^//ra  of  the  axis  from  a  straight  line; 

e'^dr  (y<l)  (16a) 

“ymax 

and 

^ebc“ 

J 

The  term  y  in  £q.  (16a)  denotes  the  position  angle  on 
the  uniformly  bent  cylinder  and  IpcYnM^L  (see  Ref.  2  for 
geometry).  The  terms  chosen  in  the  exponent  in  £q.  ( l6b) 
relate  the  curvature  of  the  axis  to  a  parabolic  curve.  Here, 
2{L/2)B^Aa,  L^fla,  and  x  is  the  position  along  the  cyl* 
inder  projected  onto  a  straight  line  (see  Ref.  7  for  geome¬ 
try).  ^  and  A  are  parameters  indicating  the  length  and 
deflection  in  terms  of  the  cylindrical  radius.  These  equa¬ 
tions  have  the  form  of  Fresnel  integrals.  As  was  done  in 
Ref.  7,  a  change  in  the  variables  to  the  following  parame¬ 
ters  is  performed  for  Eqs.  (16a)  and  (16b),  respectively: 


( rr/2 )  s  * ,  py  and  f =  yi/n )  k ,  (17a) 

and 

(iT/2)^^  =  9,k\Ax^/0^a  and  g, s2  /fa/Tr, 
which  results  in  the  compact  formulas. 

^ebc=  1 

and 

=  L  I C(^, )  +  '^( )  ] .  H8b) 

for  Eqs.  (16a)  and  (16b),  respectively,  where  Aa  is  the 
maximum  deflection  of  the  cylinder,  ^p^|/2  is  the  radius 
of  the  first  Fresnel  zone  of  the  bent  cylinder,  and  the 
Fresnel  integrals  C  and  5  are  given  as 


The  Appendix  contains  expressions  for  C  and  S  for 
arbitrary  bends  (i.e.,  deflections  are  not  restricted  with  re¬ 
spect  to  wavelength). 

A  Small  deflections.  In  the  limit  of  deflections  of  the 
bent  cylinder  that  are  sm?ll  compared  with  a  wavelength, 
the  above  equations  reduce  to  the  straight  cylinder  case 
and  ^  2. 

c.  Large  deflections.  For  deflections  of  the  bent  cylin¬ 
der  much  greater  than  a  certain  fraction  of  a  wavelength, 
or  more  precisely,  or  %k^Al}/ffla'^  1,  one  can 

use  the  method  of  stationary  phase  or  asymptotic  limits  of 
the  Fresnel  integrals  (the  Appendix)  to  solve  for  in 

Eqs.  (16a)  and  (16b); 

(20a) 

and 

Lg^^KL/4)  >jAi/Aae‘”^*,  ( 20b ) 

respectively.  Substituting  these  expressions  into  Eq.  (15b) 
gives 

/bc=*  ^/p,2,/2e'^Ac(0-O”)/^'.  (21a) 

A-i  ylA\/Aae'<’^*AJ9=01,  (21b) 

where  the  scattering  amplitude  is  shown  to  be  equal  to 
within  a  phase  factor  to  the  product  of  the  radius  of  the 
first  Fresnel  zone  of  the  bent  cylinder  and  the  normalized 
scattering  amplitude  of  the  straight  cylinder  at  normal  in¬ 
cidence  in  Eq.  (21a)  and  proportional  to  the  product  of  the 
square  root  of  the  ratio  of  the  wavelength  to  maximum 
deflection  of  the  bent  cylinder  and  straight  cylinder  scat¬ 
tering  amplitude  in  Eq.  (21b). 

II.  NUMERICAL  SIMULATIONS 

In  this  section,  the  ray  solution  is  explored  under  a 
variety  of  conditions.  The  material  properties  and  shape 


(17b) 

(18a) 
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riG  2.  Comparison  between  simple  approximate  ray  solution  [Eq.  (5)] 
illustrated  in  Fig.  1,  and  the  exact  modal-series  solution  [Eq.  (2),  math¬ 
ematically  converged]  for  weakly  scattering  infinitely  long  smooth  fluid 
cylinder  {k=k^).  Equation  (5)  is  evaluated  for  both  —tt/2  (no 

phase  correction)  and/i^^,=  —  (7r/2)k|a/(k|a-f-0.4)  (phase  correction) 
where  in  the  latter  case,  the  phase  shift  due  to  the  caustic  is  heuristically 
removed  in  the  low  ka  region.  The  ray  solution  is  remarkably  similar  to 
the  exact  solution  for  small  contrast  in  material  properties 
[(g,A)  S  (1.01,1.01),  not  shown],  especially  once  the  caustic-related  phase 
shift  is  removed  for  small  ka.  For  the  zooplankton  material  properties 
used  in  the  example  shown,  the  accuracy  is  quite  satisfactory.  Similar 
accuracies  were  observed  for  ka>  10. 


are  varied  and  compared  with  the  modal-series-based  so¬ 
lution  to  illustrate  the  usefulness  and  range  of  validity  of 
the  approximate  solution.  All  calculations  involve  back- 
scattering. 

The  weak  scatterer  approximation  is  tested  in  Fig.  2 
where  the  ray  solution  is  compared  with  the  exact  modal 
series  solution  for  values  of  speed  of  sound  and  density 
similar  to  those  observed  for  marine  organisms.  The  ray 
solution  is  given  for  cases  with  and  without  the  phase  cor¬ 
rection.  There  are  significant  difierences  between  the  ray 
solution  without  the  correction  and  the  exact  solution 
while,  once  corrected,  the  ray  solution  follows  the  exact 
solution  to  a  fair  degree  of  accuracy.  For  values  of  g  and  h 
greater  than  1.1,  the  ray  model  no  longer  is  valid  as  the 
peaks  in  the  curve  will  deviate  by  more  than  10%.  Similar 
to  the  analysis  of  Clay,^  the  impulse  response'’’^’  of  the 
form  function  of  the  exact  solution  from  Fig.  2  is  shown  in 
Fig.  3  to  decompose  the  echo  into  the  different  arrivals 
caused  by  an  incident  impulsive  (delta  function)  signal. 
The  time  of  arrival  of  the  two  echoes  observed  correspond 
to  reflections  from  the  front  and  back  interfaces  of  the 
cylinder.  The  ray  approach  is  incorporated  into  the  de¬ 
formed  cylinder  formulation  in  Fig.  4  to  estimate  the  scat¬ 
tering  by  finite  length  bodies.  Both  the  modal-series-based- 
and  ray-based  solutions  are  illustrated  for  comparison. 
Note  that  although  the  modal  series  solution  is  exact  for 
the  smooth  infinitely  long  cylinder,  the  deformed  cylinder 
solution  is  approximate,  hence  both  sets  of  curves  shown  in 
Fig.  4  involve  different  degrees  of  approximation.  Details 
of  the  calculations  involved  in  Fig.  4  are  discussed  in 
Ref.  2. 
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FIG.  3.  Impulse  response  of  exact  (modal  series)  solution  to  infinitely 
long  weakly  scattering  fluid  cylinder  (solution  taken  from  Fig.  2  with  the 
high  ka  values  tapered  via  i^e  edge  of  a  Hanning  window).  The  first 
(/p=o)  and  second  arrivals  are  illustrated. 


III.  COMPARISON  WITH  LABORATORY  DATA 

In  a  recent  study  published  by  Chu  et  al. ,  the  back- 
scatter  data  from  two  live  decapod  shrimp  were  presented 
and  analyzed.*  The  animals  were  tethered  by  human  hair 
and  individually  insonified  over  the  frequency  range  300- 
650  kHz  at  near-normal  incidence.  In  order  to  mount  the 
animals,  they  were  anesthetized,  although  note  that  they 
changed  shape  and  orientation  throughout  the  experiment 
(i.e.,  the  medication  slowed  them  down,  but  did  not  render 
them  lifeless).  Because  of  the  motion  of  the  animals,  the 
echo  level  and  spectral  shape  subsequently  changed  in 
time.  To  study  this  effect,  many  echoes  were  recorded  from 
each  animal  representing  a  range  of  realizations  of  shape 
and  orientation.  The  resultant  plots  of  backscatter  versus 
ka  indicated  strong  modal  resonances  with  nulls  ranging  in 
level  from  10  to  30  dB  below  the  surrounding  value.  The 
position  of  the  nulls  tended  to  shift  shghtly  from  realiza¬ 
tion  to  realization  as  one  would  expect  for  a  body  that  is 
undergoing  changes  in  shape  and  orientation.  Once  all  of 
the  echoes  were  averaged  (incoherently)  in  the  computer, 
the  nulls  were  mostly  “washed  out”  leaving  more  gradual 
dips  of  the  order  5-8  dB  below  local  levels. 

One  plot  representative  of  the  single  echoes  or  realiza¬ 
tions  was  compared  with  a  single  realization  of  the  modal- 
series-based  bent  cylinder  solution.  There  was  reasonable 
qualitative  comparison  between  the  structure  of  the  data 
and  theory  in  that  the  general  {ka)  position  of  the  peaks 
and  nulls  generally  agreed.  There  was  no  basis  for  compar¬ 
ison  between  the  absolute  levels  as  the  precise  shape  and 
orientation  of  the  animal  were  not  known. 

The  modal-series-based  solution  is  replaced  in  Fig.  5 
with  the  simpler,  ray-based  model.  A  similar  degree  of 
qualitative  correlation  between  the  ray  model  and  data  ex¬ 
ists  as  with  the  modal  solution  (see  above  discussion).  As 
in  the  analysis  involving  the  modal-series-based  solution,  it 
is  important  to  point  out  that  the  comparison  with  the 
ray-based  model  should  be  qualitative  also  since  the  precise 
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FIG.  4.  Ray  and  modal  series  solutions,  each  incorporated  into  deformed 
cylinder  solution,  to  estimate  scattering  by  various  deformed  elongated 
finite  length  smooth  bodies.  The  following  (zooplanktonlike)  parameters 
were  used  in  the  calculations;  (g,A)  =  (1.03S7,  1.0279)  (all  plots),  and 
Pc/L=  1.5  and  L/a=  10.5  (bent  cylinder).  The  prolate  spheroid  calcula¬ 
tions  are  valid  for  L/2ai5.  The  modal  series  is  mathematically  con¬ 
verged  in  each  plot.  Equations  (13),  (14),  and  (21a)  were  used  in  the 
(phase  corrected  or  “modified”)  ray  calculations  while  Eqs.  (7)  and 
(12b)  were  used  for  the  modal-series  curves.  Similar  accuracies  were 
observed  for  ka>  10. 


shape  and  orientation  of  the  animals  were  not  known.  It  is 
only  the  trend  and  general  positions  of  peaks  and  nulls  of 
the  data  and  model  that  should  be  compared,  not  absolute 
level  and  precise  positions. 


FIG.  S.  Comparison  between  shrimp  backscatter  data  (near  normal  in¬ 
cidence)  collected  by  Chu  et  aO  and  (phase  corrected)  ray-based  bent 
cylinder  model  [Eq.  (21a)].  Same  parameters  are  used  in  this  model  as 
those  used  in  modal-series-based  modeling  in  Chu  et  aL  (,p/L=  1.5,  L/a 
=  10.5,  and  g,)i=  1.0357,  1.0279).’  Here,  a  is  the  average  cylindrical  ra¬ 
dius  of  the  thorax  section  of  each  shrimp.  The  simple  model  predicts  the 
general  locations  of  nulls  as  well  as  overall  levels  although  note  that  the 
change  of  shape  and  orientation  of  the  live  moving  animals  caused  the 
locations  and  levels  to  shift  from  echo  to  echo  (not  shown). 


IV.  DISCUSSION  AND  CONCLUSIONS 

Perhaps  one  of  the  most  notable  results  of  this  article 
is  the  simplicity  of  the  (approximate)  two-ray  solution  and 
the  close  agreement  between  it  and  the  exact  modal  series 
solution.  Equation  (5)  shows  that  the  ray  solution  for 
weakly  scattering  infinitely  long  straight  cylinders  is  the 
summation  of  two  simple  terms  which  is  in  sharp  contrast 
with  the  much  more  complex  (exact)  modal  series  solution 
given  in  Eq.  (2).  The  simulations  show  that  for  weakly 
scattering  fluid  cylinders,  the  ray  and  series  solutions  are 
quite  close  to  each  other.  In  fact,  as  g  and  h  are  reduced  to 
near  unity,  the  ray  solution  becomes  indistinguishable  from 
the  modal  series  solution  over  a  wide  range  of  ka  (not 
shown).  In  addition,  the  real  and  imaginary  components  of 
each  solution  (not  shown)  compare  very  well  with  the 
corresponding  component  of  the  other  solution.  Once  the 
density  and  speed  of  sound  contrasts  are  increased  the  ray 
and  series  solutions  begin  to  diverge  [a  discrepancy  that  is 
not  corrected  by  removing  the  weak  scatterer  approxima¬ 
tion  in  Up  in  the  divergence  factor  in  Eq.  (3)].  Certainly 
other  phenomena  such  as  including  more  rays  would  im¬ 
prove  the  accuracy  of  those  results.^’  The  fact  that  the 
scattering  by  weakly  scattering  cylinders  was  dominated  by 
the  rays  reflected  from  the  front  and  back  interfaces  was 
illustrated  in  Fig.  3  when  the  impulse  response  of  the 
modal  series  solution  was  calculated.  The  resultant  echo 
was  composed  of  two  separate  signals  whose  times  of  ar¬ 
rival  were  consistent  with  what  one  would  expect  for  the 
two  rays.  These  times  compared  exactly  with  those  of 
Clay’  who  performed  a  similar  calculation.  When  the  ray 
solution  was  incorporated  into  the  deformed  cylinder  so- 
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lution  for  the  deformed  objects  in  Fig.  4,  the  ray-based 
solution  compared  well  with  the  modal-series-based  solu¬ 
tion. 

Another  interesting  feature  of  the  plots  is  the  fact  that, 
although  ray  solutions  typically  are  accurate  only  for  very 
high  ka,  the  "phase  corrected”  one  illustrated  in  Fig.  2 
appears  to  maintain  reasonable  accuracy  down  to  about 
/:a=0.1.  In  particular,  the  phase-corrected  ray-based  de¬ 
formed  cylinder  solution  was  accurate  down  to  values  of 
Aa=0.1,  0.02,  and  0.3  for  the  straight  finite  cylinder,  bent 
cylinder,  and  prolate  spheroid,  respectively.  As  discussed 
in  Sec.  I  A,  the  correction  involves  heuristically  removing 
the  phase  shift  induced  by  a  (high-frequency  region)  caus¬ 
tic  in  the  low-frequency  region.  Accuracy  in  this  low- 
frequency  region  is  essential  in  using  sonars  to  classify 
zooplankton. 

The  usefulness  of  the  ray-based  deformed  cylinder  so¬ 
lution  was  illustrated  in  Fig.  S  where  it  was  compared  with 
the  data  set  from  Chu  et  aL  The  simple  solution  compared 
reasonably  well  on  a  qualitative  basis  with  the  trend  and 
structure  of  the  scatter  pattern.  Certainly  it  compared  to 
the  same  degree  as  the  approximate  modal-series-based  so¬ 
lution  did  in  Ref.  5.  While  this  plot  involved  only  one 
realization  of  animal  shape,  the  data  from  other  echoes 
(some  of  which  are  shown  in  Ref.  S)  contained  similar 
structure.  In  particular,  all  data  contained  peaks  and  nulls 
(or  dips)  occurring  in  the  same  general  location  as  that 
predicted  by  the  single  realization  predictions.  The  shifting 
of  their  locations  from  echo  to  echo  is  undoubtedly  due  to 
the  constantly  changing  shape  and  orientation  of  the  ani¬ 
mals.  A  more  rigorous  analysis  in  predicting  the  precise 
locations  of  the  peaks  and  nulls  and  absolute  levels  would 
require  knowledge  of  the  precise  shape  and  orientation  of 
the  animals.  Nonetheless,  the  favorable  comparison  sug¬ 
gests  that  the  animals  (decapod  shrimp)  can  be  treated  as 
weakly  scattering  fluid  bodies  whose  scattering  properties 
near-normal  incidence  can  be  approximated  by  a  simple 
two-ray  model. 

In  addition  to  the  major  peaks  and  nulls  mentioned 
above,  there  also  appear  to  be  some  smaller  ones  occurring 
near  ka=2.3,  3.2,  4.2  (possibly),  and  5.4.  Since  the  struc¬ 
ture  is  so  weak  and  the  data  set  so  limited,  explanation 
regarding  the  mechanism  causing  the  oscillations  would  be 
purely  speculative.  However,  their  near-unity  periodicity 
does,  in  fact,  correspond  to  the  periodicity  one  would  ex¬ 
pect  due  to  interferences  between  circumferential  waves 
traveling  near  the  speed  of  sound  of  the  zooplankton  and 
surrounding  water.  Whether  of  not  such  waves  exist  with 
these  animals  remains  to  be  seen. 

In  conclusion,  the  simple  two-ray-based  deformed- 
cylinder  solution  has  proven  to  be  a  good  descriptor  of  the 
scattering  by  idealized  weakly  scattering  fluid  objects  near¬ 
normal  incidence.  The  model,  besides  having  the  advan¬ 
tage  of  illustrating  the  basic  physics  of  the  scattering,  is  in 
such  a  form  that  can  be  manipulated  analytically  and  is 
easy  to  compute.  The  reasonable  qualitative  comparison  of 
the  structure  of  the  ray  solution  with  the  shrimp  data  sug 
gests  that  the  irregularly  shaped  animals  can  be  modeled  as 
weak  fluid  scatterers  whose  scattering  properties  can  be 


approximated,  at  least  near-normal  incidence,  by  a  simple 
two-ray  model.  Although  the  scattering  structure  pre¬ 
dicted  by  the  single  realization  ray  model  is  very  sensitive 
to  the  precise  shape  and  orientation  of  the  animals,  the 
sensitivity  is  dramatically  reduced  once  ensemble  averages 
with  respect  to  length  and  orientation  of  the  extended  so¬ 
lution  are  performed  in  Ref.  1 1 .  In  that  companion  article, 
the  deep  nulls  from  the  single  realization  solution  become 
rounded  “dips”  due  to  the  averaging.  The  averages  com¬ 
pare  very  well  with  data  involving  aggregations  of  up  to 
hundreds  of  animals.  Hence,  the  end  result  in  that  analysis 
involving  many  animals  further  validates  use  of  the  ray- 
based  formulas  to  describe  scattering  in  the  single  realiza¬ 
tion  case. 
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Numerical  expressions  for  the  evaluation  of  Fresnel 
integrals  are  in  Abramowitz  and  Stegun,  Sec.  7.3.^*  We 
repeat  the  expressions  here  for  the  readers’  convenience. 
The  notation  for  the  real  variable  x  is  not  related  to  x  used 
in  the  main  text.  The  equation  numbers  A&S  are  those  in 
the  handbook. 

C(x)=  J*  cos^t^^dt  (A&S  7.3.1) 

and 


4  V/. 

(A&S  7.3.32) 


APPENDIX:  CONVENIENT  FORMULAS  FOR  FRESNEL 
INTEGRALS 
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5(x)=  sin^^r^jrfr.  (A&S  7.3.2) 

For  real  variables  and  using  the  auxiliary  functions 
fix)  and  g(jc),  the  Fresnel  integrals  are 

Cix)  =  2+/(-*)sin|^  j  -g(x)cos^“  j. 


(A&S  7.3.9) 


1 


5(x)---/(x)cos^-x^j-g(x)sin|"X^j. 

(A&S  7.3.10) 

Rational  approximations  for  the  auxiliary  functions 
(0<x<  00 )  are 
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'  ^  ^  1 

*^-*^=2  +  4.142x  +  3.492x^  +  6.670cr^'*'^^'*^’ 

|€U)|<0.002.  (A&S  7.3.33) 
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